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Abstract 

A curve of applied pressure P versus lattice constant a is calculated 
for single crystal aluminum. It results from an application of the method of 
structural expansions for deriving the energies of simple metals, a method 
known to give reasonable results for the elastic constants even at second 
order in the effective electron-ion interaction. The latter (in the present 
calculation) is taken from Fermi surface analysis and it is verified (with 
this essentially experimental information) that the extant face centered cubic 
structure remains the preferred crystalline phase up to the highest pressures 
considered. Arguments are given to suggest that the P vs. a curve should 
have reasonable a prior i accuracy, and can admit of possible improvement if 
experimental data in the intermediate pressure region can be provided to 
refine the (in principle) energy dependent pseudo-potential. At 3 megabars 
the lattice constant is reduced by only 22%; the ion cores at this pressure 
are still very well separated. 
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l. 1 n tvodttc t ion 

Among the simple metals, aluminum is in many ways one oT the simplest 
being cubic close packed under normal conditions and possessing, ion cores 
occupied by electrons in levels of s- and p- symmetry . It is Plainly a conse- 
quence of the latter that its nearly f ree-electron band structure can be 
interpolated so accurately by a spatially local pseuilopotential , a feature 
which distinguishes it somewhat from the alkali metals. Although the Fermi 
surfaces of the alkali metals are a good deal simpler than that of aluminum, 
the apparant complexity of its multiply-connected Fermi surface can be used 
to advantage in a study of the transport properties at high pressure. This 
will be the content of a later work; for the present we are concerned with 
the equation of state of Al, a necessary preliminary in discussing the depen- 
dence of transport properties on pressure. 1 Effects of temperature (for 
normal conditions) are quite small and our nim here is therefore to express 
the equation of state in terms of pressure versus lattice constant. Such a 

relation can only be considered potentially useful If no crystalline phase 

2 

changes are likely to occur. We show by a series of arguments that the 
common face centered cubic phase of Al appears to remain the stable phase for 
pressures exceeding 3 megabars. In terms of the lattice constant (or equi- 
valently the r g electron spacing parameter) these collosal pressures 
represent a rather modest change of around 20%. The electron density is 
increased, but not greatly. It is not unreasonable to suppose therefore that 
the method based on structural expansions about the uniform interacting elec- 
tron gas will continue to function ns it does for the system taken at more 

reasonable pressures. The method is summarized in the noxt Ruction (II), and 

3 

in the course of discussing the standard second order theory we comment, on 
the importance of higher order corrections to the present calculations. 





Section III describes- the application of the formalism to Uie problem of 
deciding which of several possible simple structures (including fee) will 
possess the lowest Gibbs energy. For the fee phase, a curve of lattice 


constant a versus pressure is presented (Section IV); up to and above 3 nhar, 

the changes in a nve quite mono tonic. We estimate. that not until pressures 

of over 100 mbar are reached will the ion cores of.Al bo substantially 

contiguous. This is a very different situation front the one prevailing in 

ionic crystals where the pressure scale is founded largely on assumed short 
4 

range interactions. Although the atomic number of A1 is relatively low it 
may compete reasonably well in X-ray scattering power with NaCl and may. 
therefore be an alternative candidate for calibration and use as a pressure 
scale. 


t- ... 
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H. Energy of Simple .Metals 


On account of the compactures of its ion cone (anil the absence of 

filled d-shell levels) the pseudopotent.ini in Al, although energy dependent to 
(5) 

a small degree is remarkably local and provides an excellent interpolation 

to a priori bend structures. Invoking an adiabatic approximation, fra shall taka 

it that an ion of the dynamic lattice of Al carries with it a bare pseudopotential, 

(C s 

v(k) known (at the Fermi energy) from Fermi surface analysis ; It is a 

function which as is well known oscillates in sign as k increases, a fact 

which reflects the finite size of the Al ion cone. Since we shall shortly 

need to consider the possibility of corrections arising from dynamic lattice 

effects, it is convenient to set down a Hamiltonian for the electron system 
(7) 

that is written for Instantaneous portions r(R) of the ions near equilibrium 


sites R, l.e. 


II a H + II + H 
eg M ei 


where for the present H can be taken as the standard Hamiltonian for the 

eg 

interacting electron gas (uniform compensating positive background) and the 

ionic Hamililtonian II leads to the Madelung energy NE of point ions. In 

M M 

Rydbergs it can be written (for ZN electrons) 

:Y 2NE H =i&nr <*<**•» >> 

where for the ion3 in a volume V the structure factor for the ionic system is 


s <« - 5 'W' “k.o 


\ - 5 


,ik • r (R) 


and the average in (3) being taken over the states of the crystal. The final 
term in (1), H , is the electron ion interaction in which it is convenient to 
include the largely compensating zeroth Fourier component of all the long range 
interactions', that, is, a term which although independent of structure is 
Always difficult to calculate from first principles. It can, however, be 
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el irai noted by exploiting :i f ragmeut of experimental information such as, for 
example, the equilibrium density . ^ ^ 

Accordingly we write 


It . = K +.5’ ft v(k) '? ® 
ei n kgO k — -k 

where for the electrons the density operator is written 


ac . , 

(i = L e 

- i 


ik • r 


( 1 ) 


C5) 


We turn first to the static lattice case for which the contribution of 

E (] to the thermodynamic functionals known, at least for most simple structures. 

The problem of calculating the energy of a simple metal then reduces to an 

expansion (relative to the structueless electron gas system) in orders 

(beginning at the second) of J| . Since the ionic correlation function (for 

example, S(k)) are then de.lta functions of the reciprocal lattice they reduce 

the resulting summatious in the perturbation series to lattice sums. Thus, 

(9) 

in addition to the ground state energy from II . (and E, t ) we have, as the 

eg M 

first term the structural expansion a second order contribution of the form 
E = l y' |v(k)j 2 (k)/e(k) ; [k] reciprocal lnttice set (6) 

EQ no 


where e (k) is the dielectric function of the interacting electron gas and 
X a) (k) its (static) first order polarizabitity . At this level of approximation 
the internal energy is then: 

E - [E F«> <” 

eg M o 1.5 

and it i3 interesting, before proceeding further to examine their relative 

contributions to the pressure at a given volume V, or what is equivalent, a 

mean electron spacing r g ( V/NZ n (r^ao) 3 dn/3). Table X shows that as 

pressure increases the contribution from E becomes progressively a smaller 

JJS 

fraction of the total. Since we kuow^^ the ground state energy and compresibt 1 ity 
of A1 to be quite well given near P = 0 by (7) and its derivatives, we may 
conclude that even at high pressures the higher order hand structure contrtbu— 
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I'.ioiis l:o )•: are not likely to lx; :m important factor in limiting tho :ici:ur:tcy 

of a calculation of I’ vs a. Tho n;i>sL significant of those; ruiacctiims 

the third order band structure; energy, if the electron gas is treated, for 

711 

example, v.ithin the random phase approximation this term can ho written ’ 

fl/3U v(k) v(lQ v (l:-U * ) x (2) 0:, k , ) K-K * ) (8) 

kk* F(U) E (k 1 ) e (k-k ' ) 

(2) 

ivliore is the second order polarizability of the electron system. As 

remarked earlier, VOOfor A1 (and -indeed any non-point ion-system) alternates 

in 6ign as its argument increases and as a consequence there is substantial 

self cancellation in (8). Furthermore, relative to E^., the V(k) are 

considerably less than ~0.1 (for example |v, , , /E„ I =0.020!), and I V„_.,./E„|-= 

1 111 r • 1 200 F . 

0.0657). It follows that the higher order band structure energies nre quite 

(2) (12) 
small in comparison with E^ , This lias already been noticed by others . , 

although we must ■ recognise that the derivatives of the higher order terms (in 
the clastic constants for example) need not always be unimportant. 

As far ns a calculation of the pressure is concerned it seems a reason- 
able approximation to neglect the higher order hand structure energies. The 
approximation would appear less Justifiable in the calculation of the ground 
state energy for various crystal structures. But in fact it remains numeri- 
cally valid. The concern is that differences in Gibbs energy for different 
crystal structures are quite small, about 4-G mtly between hep. and fee per 
electron if calculated with a second order expression. And these can be less 
than typical third order energies. However, . v;e need not tho a hsolu to third 
order energies , but their differences for different structures; these arc in 
turn smaller by about an order of magnitude. Vic shall see in n moment that 

inclusion of dynamic effects. are likely to reduce the third order differences 

' i . ■ ■ . v ' 

still further, so that a calculation of the energy at second order is sufficient 
for. the present purposes. . - .1 . 



/ : / 





Hola-xing the staLic lattice assumption requires (a) llio inclusion ol 
phonon energy tens, IT indeed tlu; excitations arc to be described by phonons, 
and (b) the ' reintroduc lion in (R) and (7) of the corresponding ionic correla- 
tion functions, for example S(k) (equation (3)). If u(K) is the displacement 
Of hn ion from site K then 



1C „ e«s. H e i a .n* e - i( a,+W .JB". 
nu'u 

exp -l <(k. u(n)) 2 + (q. U ( R ’)) 2 + ((qlk) . u(R)")) 2 
j ; + 2(k ► u(R) q . u(R*)) - 2(a . u(n') (k+q) . u(U”)) 

^ AO ro AWA ^ A/ AW* AWA # A AXA 

r;- -2(k.u(ft) (k+q) . u(R"))> 

^AWAW AW *wA . ^ , • 

.which is straightforward to generalize to higher orders. 

For. metals with substantial Debye temperatures (in which category we may 
’ place Al) one method of handling (10) and (11)) is to proceed by a multi- 
phonon expansion. Tlie zero phonon term leads immediately back to (6) and (8). 

. The one-phonon term leads, when combined with the kinetic energy of the phonon 
system ^ to the internal energy of the phonons. The remaining multi plionon terms, 
.-'us is known from the analysis of thermal diffuse x-ray scattering arc quite 
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small. Thus wo may, with n sufficient accuracy. Iron t the phonons independently 
of the electron system and calculate the Gibbs energy of the latter a.ssumiug 
a rigid lattice. The internal energy can then be written 


E = (v: + E + K 1 + F,j 2) l- 

V eg o M ./ BS 

where is the internal energy of the phonon system. 


( 12 ) 



I 







III. Structural Consi derat ions 

From the known l-’orm l surface of A1 (ami the assumption of ;i static laLl'ico) 

the values of v(X), K-(lll), (200) can he extracted and these can he Inter- ■ 

o 

polated and extended by an empty-core pseudopotential (v(K) = C-fcnZ/jOcos |lcjr ) 
The rnnge of validity (in k) of such a simple. form is quite sufficient to 
assure convergence of the sums in (0), and hence of the band structure energy. 
Since v(K) is a property of the ion we nay repeat the procedure at any chosen 
volume or density. Assuming for the moment that this is fixed we must examine 
the structure dependent terms in (12) as the ions are rearranged in a variety 
of possible crystal structures. 

To begin with we consider. the electronic terms (and Madelung energy) ard 
allow ourselves at this point the freedom of a structure with a two atom basis. 
Tiie task is to ascertain which of. the structures (at least, which of the simple 
structures) is preferred for Al: to this end we will select carefully n. system 
of primitive and basi3 vectors which will allow us continually to deform 

15 

between different structures by means of a smooth variation of parameters. 

Refer now to Figure 1(a). We take a, b and c ns primitive vectors which are 
written in the form: 

a = a(s,o,o) , 

b - n(v* ,f;,o), (13) 

c = a(o,o,ri) . 

. *** 

Direct lattice vectors are then written 

R » na + pb + qc . . 

We take the basis vectors 

• *1 = 0 » ' 

~2 * I = f (2s " 1 * <2s -Dc» t i ) - do 

In (13) and (14) the parameters v’, s, q, and £ arc chosen in the followiiR way 



1(1 - 


v' (2s-l ) v 

§ = u-v(2s - /3) i (1-s) l’l-2u ^ 2vl2u-/3}l 

H = v + 2 vw (/sTTt - 1) i 2 (1-s) U-w-2vw (/'S/Tt-l ) I 
£ = U-2r (u - \/J 3 ( 10 ) 

with s,u,v,w, taker, as i tide pend eat parameters. Transformation (15) is only 
one of many ways of continually deformin'' the standard simple crystal structures. 
We have selected it because it permits us to examine single cubic (sc) face- 
centered cubic (fee), body centered cubic (bcc), and hexagonal close packed (hop) 
with variable (c/a) ratio. As an example, note that when s = J we have (what- 
ever finite values u,v,w may assume) a simple cubic structure. Cn the other 
hand, if 3 = 1, v=G, and w = 1, the structure is fee for u sc/2, and bcc with 
u = 1 . Further, if s = 1, v = and w = 1. we have hep with ideal ratio-. ITiese 


are summarized on Figure 1 (b) and (c) . Although it cannot be deduced simply 
from the results we shall give, it is interesting to. note that the transformation 
we have chosen moves the atoms in a very natural way, keeping them well apart, 
and proceeding as directly ns possible from one structure to another. In a sense 
we are moving the atoms along valleys in the energy-structure space. 

The lattice reciprocal to (13) is spanned by primitive vectors 


& ■ “ ( b- % ■ °> • 

E ■ T <0 - V 0> . 


_ 2TT ■ 1 . 

C = — (0, 0, ) , 


and the reciprocal lattice vectors are 



K = h A + £B + tnC. 
which we use to define in A1 (Z=3) 


(2k ) _1 K 
F ~ 


* = < ~ s 5 n) 2/3r “ 2 


L-2 < 


hv* 


2 -.ft 

m | 

+ 7F J- 


(17) 


With the choice of basis given in (13) the structure factor, per ion, is 

} O ( e -ig ) . ■■■'- . 


where 




thnn 



I 


c “ V u “* 3> " 71 ' _t 


I* ( i:i -r<: ('.l-t ) cm Ic(^!>_) 


crfc ( ■ 1 .I»|.- h -L|) 
|p + t| 


where r is Kwald's dimensionless parameter and erfc denotes the complementary 

( 2 ) 

error function. Then at second older, v/e evaluate (11) by usinsc (r») for F. 

■ 

- — i fi 

(with v(l<) there replaced by 1 1- e (k)) and (■■(9) amt (20) for tlie Madelung 

2 ■ 

energy)- For a given structural choice (corresponding to a particular selection 

of s,u,v,wj we determine E 0 by the zero pressure condition (JK/Jr ) _ =0. 

s r 

. so 

Expressed as an energy per electron, E always has the form 

, , 4tt 3. 
u/( — r ) 

3 s 

where u Is a property of the ion alone and is assumed not to alter under reasonable 
variations of density. Since the total energy near 1’ = 0 actually contains 
small contributions from the omitted .higher order band structure terms, the 
imposition of the zero pressure condition forces their inclusion in a crude 
way through the choice of a. To the extent that these terms are not seriously 
density dependence the subsequent use of this a will therefore continue to 
incorporate such terms. If one takes the h'ozieres -Pines form for the correlation 


energy ; it is easy to see that 


u = ~ [r 2 (0.916 + Z 2/3 C , + 0.031 r + r 2 ) ) - 4.42r 1 (21) 

9 so M so so \<> rs / so 

■ ■ ' • r 

so 


where for the fee structure observed for A1 in its ground state r = 2.0647. 

so 

(21 

What is required in (21) is (r ) nnd this can be calculated bv a combination 

of a direct numerical summation (out to a chosen reciprocal lattice shell) 

augmented by. integration for the remainder. This remainder, designated by 

S(x x , r s ) (where x is the radius of the shell) is independent of structure. 

18 • 

and depends . very weakly on r^ . Its contribution is in any event, quite 

small. At r = r and for x, = 2.5 we find S = 0.005 ly/electron which amounts 
s so 1 

( 2 ) 

to 55' of K and 0.41- of K. 

I*I> 



IV. Energies and Phases: Results and Discuss! oils 

In figure 2 we show a selection of the results v;e obtain Lor the Ilolmol Lz 
free energy E as the crystal structure is continuously deformed from fee to 
hep (c/a = /8/3) . In this example fee is lower in energy at all densities 
considered. This result remains true for other structures, the two that are 
always closest in energy (at least of the simple structures we consider) 
being fee and hep. It is a straightforward matter to compute the PV term 
and hence, in the ground state the Gibbs energy for different phases. We find 
fee A1 (with an assumed static lattice) to have the lowest Gibbs energy and 
to be the preferred structure, even up to theoretical pressures in excess of 
3 mbars. 

Contributions to the thermodynamic functions from the ionic degrees of 

freedom can be estimated from the Debye model; in particular, the zero point 

energy is of order (3/8)k„0 r , per-electron (about 0.001 Ily) and for tempera- 

B U ' 

tures less than the Debye temperature will remain of this order. Changes in 
this energy accompanying changes in crystal structure will be much less than 
0.001 Ry. The contribution of the phonons to the pressure is readily shown 
to be where n = N/V is the ionic density, and y = the Griineisen 

constant. Even for changes of 50% in the equilibrium value of n, the phonons 
change the pressure calculation above by at most a few kilobars. Figures 3 


and 4 give the Gibbs energy as a function of pressure for fee and hep, and 
(for comparative purposes) as a function of r^ fc sc, fee, bcc. and hep. 

- In figure 5 we plot the pressure on a single crystal of A1 (under pure hydro- 
static strain) as a function of its lattice constant a (rather than r g ) at a 
nominal temperature of 300°K. The equation of state given there may also bo 


appropriate to polycrystalline samples under less that pure hydrostatic 

. • o 

conditions. It is worth remarking that at 3 mbnr, where a = 3.14 A, and the 


nearest neighbor separation is - 777 a ~ 2.22 A, the distance between ion cores 






e t ' 

(taking them to have n radius oC 0.59 A) is still 1.01 -A. For the pressure 
range- given in figure 5 the energy (and the corresponding pressure) is dominated 
by the terns arising from electron gas. Modelling energy, and to a much lesser 
extent, band structure. Energies arising from the direct overlap of ..ion-cores 
(so called core-core exchange, or Born-Mayer terms) are evidently not important 
although it is conceivable that at very much higher pressures (ive estimate they 
will be in excess of 100 Mbars) they could be. This kind of torn is difficult 
to calculate with confidence from first principles, and 'is normally para- 
meterized- in. an exponential form (or even ns a power law) in expressions giving 
its contribution to the internal energy. In pressure scales based on these 
forms, the concern (aside from the implicit pair force approximation) i3 that 
the low pressure determined parameters may not remain valid in a region of 
substantial ion-core wave function overlap. At 3 mbar we have only a 22% 
reduction in lattice constant and core-core overlap is still a small effect; 

its neglect leads to errors which will be far less important than those 
arising from the neglect of, for example, the higher order band structure 
contributions to the energy. 

As for as the use of A1 in devices is concerned it suffers from the dis- 
advantage that its atomic number is quite low. It- should, however, be visible 
to x-rays in a diamond cell and the curve presented in figure 5 is therefore 
amenable to experimental test, provided of course that sufficiently hydrostatic 


conditions can be arranged. ' ■ • -;r j; .V .-...vv 

If a test of this kind were found to establish ns numerically sound the 

- ‘ " . ' ' o'. » ' 

bas 4 c . curve up to, say, 0.5 m bar (corresponding to a = 3.01 A) then according 

to the agruments wo have given above it would then appear reasonable to accept 

19 

the balance of the curve leading to the ultra pressure region. An independ- 

ent determination of the pressure can also be used to refine, for example, the 
form of the pseudopotuntial used in the high density regime. Finally the 

reasons fin- choosing A1 (the paradigm. of small core, close-packed— cubic nearly 
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/ Figure Captions 

Figure 1: (a) General structure defined. (b) Some particular cases and represen- 
tations of continuous one-parameter transformations of them into each other. 

(c) Values of the parameters for these particular cases. The parameters nre • 
defined by Eqs. (13) - (15). 

Figure 2: Helmholtz free energy as a function of and v, the other para- 

meters fixed at their fee values; varying v here takes the structure from fee 
to hep. ; 

Figure 3: Gibbs free energy as a function of r s for several common structures. • 

Compare with Figure 4 where G is plotted against the natural variable, pressure. 

Figure 4: Gibbs free energy a3 a function of pressure for the fee and hep 

structures; these have the lowest Gibbs free energy for any fixed pressure P. 

Figure 5; Pressure as a function of lattice constant for the fee structure. i 

2 3 

An approximate fit for it is P = 0.786 (a^-a) - 0.824 (n^-n) +4.131 (a o ~a)‘ , 

the maximum error being 13% at about 0.1 Mbar, but less that 4% for most of I 

the range of a. | 






Table ( 


r 

P[K 

+ K + E 1 

(2) 

l‘[K~] 

s 

eg 

M o 

ns 

2.07 

O.'JS 

[-1.29] 

-0.18 [-0.097] 

1.9 

1 .39 

[-1.24] 

-1.07 [-0.133] 

1 .8 

2.38 

[-1 .189] 

-1.62 [-0.176] 

1 .7 

3.95 

[-1.110] 

-2.37 [-0.227] 

1.6 

6.47 

[-0 . 993 ] 

-3.35 [-0.292] 


Caption: The quantities E , K and E arc present at any order of the 

eg M o 

calculation and are convenient to group together in the comparison of the 

relative pressure contributions. The first column gives an estimate of the 

( 2 ) 

pressure (in Mbars) from [E + E +E ] and the second column for E__ . 

eg M o ES 

Energies are given in Rydbergs. 





































